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Abstract 

We report on experiments in 6 Li Fermi gases near Feshbach reso- 
nances. A broad s-wave resonance is used to form a Bose-Einstein con- 
densate of weakly bound 6 Li2 molecules in a crossed optical trap. The 
measured molecule-molecule scattering length of 170lgQ° nm at 770 G is 
found in good agreement with theory. The expansion energy of the cloud 
in the BEC-BCS crossover region is measured. Finally we discuss the 
properties of p-wave Feshbach resonances observed near 200 Gauss and 
new s-wave resonances in the heteronuclear 6 Li- 7 Li mixture. 

Strongly interacting fermionic systems occur in a variety of physical processes, 
ranging from nuclear physics, to high temperature superconductivity, superfluidity, 
quark-gluon plasmas, and ultra-cold dilute gases. Thanks to the phenomenon of 
Feshbach resonances, these gases offer the unique possibility to tune the strength and 
the sign of the effective interaction between particles. In this way, it is possible to 
study the crossover between situations governed by Bose-Einstein and Fermi-Dirac 
statistics. 

1 BEC-BCS crossover near 6 Li s-wave resonance 

When the scattering length a characterizing the 2-body interaction at low tempera- 
ture is positive, the atoms can pair in a weakly bound molecular state. When the 
temperature is low enough, these bosonic dimers can form a Bose-Einstein conden- 
sate (BEC) as observed very recently in 40 K [I] and 6 Li EI EI El- On the side of the 
resonance where a is negative, one expects the well known Bardeen-Cooper-Schrieffcr 
(BCS) model for superconductivity to be valid. However, this simple picture of a 
BEC phase on one side of the resonance and a BCS phase on the other is valid only 
for small atom density n. When n\a\ 3 > 1 the system enters a strongly interacting 
regime that represents a challenge for many-body theories 0. In the recent months, 
this regime has been the subject of intense experimental activity |4ll^ll7ll51l^l llUltTT| . 

Here we first report on Bose-Einstein condensation of 6 Li dimers in a crossed op- 
tical dipole trap, and a study of the BEC-BCS crossover region. Unlike all previous 
observations of molecular BEC made in single beam dipole traps with very elon- 
gated geometries, our condensates are formed in nearly isotropic strongly confining 



traps. The experimental setup has been described previously |12l 1 1 3 1 . A gas of 
6 Li atoms is prepared in the absolute ground state |l/2, 1/2} in a Nd-YAG crossed 
beam optical dipole trap. The horizontal beam (resp. vertical) propagates along x 
(y), has a maximum power of P„ = 2W {P„ = 3.3 W) and a waist of ~ 25 /tm 
(~ 40 /im). At full power, the 6 Li trap oscillation frequencies are lu x /2tt — 2.4(2) kHz, 
LO y /2n = 5.0(3) kHz, and u z /2n = 5.5(4) kHz, as measured by parametric excita- 
tion, and the trap depth is ~ 80 (mK. After sweeping the magnetic field B from 5 G 
to 1060 G, we drive the Zeeman transition between |l/2, 1/2} and |l/2, — 1/2} with 
a 76 MHz RF field to prepare a balanced mixture of the two states. As measured 
recently |11|. the Feshbach resonance between these two states is peaked at 834(2) G, 
and for £>=1060G, a — — 167 nm. After 100 ms the coherence between the two states 
is lost and plain evaporation provides JVf = N± = N to t/2 — 1.5 x 10 atoms at 
10/iK=0.8T F , where k B T F = h 2 kl/2m = h{2>N tot Lu x u y Lu z ) 1/;i = to(37V tot ) 1/3 is the 
Fermi energy. Lowering the intensity of the trapping laser to 0.1 Po, the Fermi gas is 
evaporatively cooled to temperatures T at or below 0.2 Tf and iVtot ~ 7 x 10 4 . 

Then, sweeping the magnetic field to 770 G in 200 ms, the Feshbach resonance is 
slowly crossed. In this process atoms are adiabatically and reversibly transformed 
into cold molecules 1131 I I 1| near the BEC critical temperature as presented in figure 
la. The onset of condensation is revealed by bimodal and anisotropic momentum 
distributions in time of flight expansions of the molecular gas. These images are 
recorded as follows. At a fixed magnetic field, the optical trap is first switched off. 
The cloud expands typically for 1 ms and then the magnetic field is increased by 
100 G in 50 fis. This converts the molecules back into free atoms above resonance 
without releasing their binding energy [3]. Switching the field abruptly off in 10 fis, 
we detect free 6 Li atoms by light absorption near the D2 line. We have checked that, 
in the trap before expansion, there are no unpaired atoms. In figure lb, a Bose- 
Einstein condensate of 7 Li atoms produced in the same optical trap is presented. The 
comparison between the condensate sizes after expansion dramatically reveals that the 
mean field interaction and scattering length are much larger for 6 Li2 dimers (Fig. la) 
than for 7 Li atoms (Fig. lb). 

To measure the molecule-molecule scattering length, we produce pure molecular 
condensates by taking advantage of our crossed dipole trap. We recompress the hori- 
zontal beam to full power while keeping the vertical beam at the low power of 0.035 Po 
corresponding to a trap depth for molecules U = 5.6 /iK. Temperature is then lim- 
ited to T < 0.9 /xK assuming a conservative r/ — U/ksT = 6, whereas the critical 
temperature increases with the mean oscillation frequency. Consequently, with an 
axial (resp. radial) trap frequency of 440 Hz (resp. 5kHz), we obtain T/Tq < 0.3, 
where Tq = to(0.82iVtot/2) 1/3 =2.7 (jK is the non interacting BEC critical tem- 
perature. Thus, the condensate should be pure as confirmed by our images. After 
1.2 ms of expansion, the radius of the condensate in the x (resp. y) direction is 
R x = 51 /im (R y = 103 /im). The resulting anisotropy R y /R x — 2.0(1) is consis- 
tent with the value 1.98 EH predicted the scaling equations Moreover, this 
set of equation leads to an in-trap radius R x = 26/im (resp. R y — 2.75/im). We 
then deduce the molecule-molecule scattering length from the Thomas-Fermi formula 
Rx, y = aho<i'/^a:,j;(15A'totam/2aho) 1 ' /5 , with aho = \/ h/2muj. Averaging over sev- 
eral images, this yields a m =170tgo°nm at 770 G. Here, the statistical uncertainty 
is negligible compared to the systematic uncertainty due to the calibration of our 
atom number. At this field, we calculate an atomic scattering length of a = 306 nm. 
Combined with the prediction a m — 0.6 a of ^S], we obtain a m — 183 nm in good 
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agreement with our measurement. For 7 Li, we obtain with the same analysis a much 
smaller scattering length of 07=0.65(10) nm at 610 G also in agreement with theory 

EH 

The condensate lifetime is typically ~300ms at 715 G (a m = 66 nm) and ~3s at 
770 G (a m = 170 nm), whereas for a = —167 nm at 1060 G, the lifetime exceeds 
30 s. On the BEC side, the molecule-molecule loss rate constant is G = 0.26to'o<3 x 
10" 13 cm 3 /s at 770 G and G = 1.75±%* x 10" 13 cm 3 /s at 715 G with the fit procedure 
for condensates described in |20|. Combining similar results for four values of the 
magnetic field ranging from 700 G to 770 G, we find G oc a ~ 19±0 - 8 (figure 2). Our 
data are in agreement with the theoretical prediction G oc a~ 2 ' 55 of ref. |18| and with 
previous measurements of G in a thermal gas at 690 G 13J or in a BEC at 764 G [B]. 
A similar power law was also found for 40 K fTTI. 

We then made an investigation of the crossover from a Bose-Einstein condensate 
to an interacting Fermi gas (Fig. Qc and d). We prepare a nearly pure condensate 
with 3.5xl0 4 molecules at 770 G and recompress the trap to frequencies of uj x — 
2-7T x 830 Hz, ujy = 2ty x 2.4 kHz, and lu z — 2tt x 2.5 kHz. The magnetic field is then 
slowly swept at a rate of 2G/ms to various values across the Feshbach resonance. 
The 2D momentum distribution after a time of flight expansion of 1.4 ms is then 
detected as previously. As B increases from the regime of weak interactions the 
condensate size gradually increases towards the width of a non interacting Fermi gas. 
Nothing particular happens on resonance. Fig. 0c and^d present respectively the 
anisotropy of the cloud after expansion r] and the corresponding released energy E Te \. 
These are calculated from gaussian fits to the density after time of flight: E le \ = 
m(2ay + a x )/2r 2 and i] = a y /a x , where m is the rms width along i, and r is the 
time of flight. The anisotropy monotonically decreases from ~1.6 on the BEC side, 
where hydrodynamic expansion predicts 1.75, to 1.1, at 1060 G, on the BCS side. 
On resonance, at zero temperature, a superfluid hydrodynamic expansion is expected 
|22| and would correspond to r\ = 1.7. We find however r\ — 1.35(5), indicating a 
partially hydrodynamic behavior that could be due to a reduced superfluid fraction. 
On the a < side, the decreasing anisotropy would indicate a further decrease of the 
superfluid fraction that could correspond to the reduction of the condensed fraction of 
fermionic atom pairs away from resonance observed in [7||5|. Interestingly, our results 
differ from that of ref.[23j where hydrodynamic expansion was observed at 910 G in a 
more elongated trap for T/Tf — 0.1. 

In the BEC-BCS crossover regime, the gas energy released after expansion E Te \ is 
also smooth (Fig. 0d). E rc \ presents a plateau for B < 750 G, and then increases 
monotonically towards that of a weakly interacting Fermi gas. The plateau is not 
reproduced by the mean field approach of a pure condensate (dashed line). This is 
a signature that the gas is not at T = 0. It can be understood with the mean field 
approach we used previously to describe the behavior of the thermal cloud. Since the 
magnetic field sweep is slow compared to the gas collision rate |13l . we assume that 
this sweep is adiabatic and conserves entropy |24| . We then adjust this entropy to 
reproduce the release energy at a particular magnetic field, B — 720 G. The resulting 
curve as a function of B (solid line in Fig.Qd) agrees well with our data in the 
range 680 G < B < 770 G, where the condensate fraction is 40%, and the temperature 
is T ~ Tq/2 = 1.4/iK. This model is limited to rimd^ < 1. Near resonance the 
calculated release energy diverges and clearly departs from the data. On the BCS 
side, the release energy of a T — ideal Fermi gas gives an upper bound for the 
data (dot-dashed curve), as expected from negative interaction energy and a very 
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cold sample. This low temperature is supported by our measurements on the BEC 
side and the assumption of entropy conservation through resonance which predicts 
T = 0.1 T F on the BCS side 

On resonance the gas is expected to reach a universal behavior, as the scattering 
length a is not a relevant parameter any more |25|. In this regime, the release energy 
scales as E rc \ = y/l + f3E® cl , where £* c i is the release energy of the non-interacting gas 
and (3 is a universal parameter. From our data at 820 G, we get /3 = —0.64(15). This 
value is larger than the Duke result j3 = —0.26 ± 0.07 at 910 G |23| . but agrees with 
that of Innsbruck (3 = — 0.68 j^Q '^q at 850 G [Sj, and with the most recent theoretical 
prediction (3 = -0.56 [101111 Around 925 G, where a = -270 nm and (k F \a\y L = 
0.35, the release energy curve displays a change of slope. This is a signature of the 
transition between the strongly and weakly interacting regimes. It is also observed 
near the same field in through in situ measurement of the trapped cloud size. 
Interestingly, the onset of resonance condensation of fermionic atom pairs observed in 
40 K [7] and 6 Li [Hj, corresponds to a similar value of /cf |a| . 

2 P-wave resonances 

Recently, both experimental 28 29 30! and theoretical |31| papers have devoted 
interest to the p-wave Feshbach resonances. The goal of these experiments is to 
nucleate molecules with internal angular momentum 1 = 1, that could lead to the 
observation of some non-conventional superconductivity, analogous to that observed 
in superfluid 3 He |55) . 

In the manifold / = 1/2 corresponding to the hyperfine ground state of 6 Li, coupled 
channels calculations have demonstrated that three Feshbach resonances could be 
observed in p-wave channels. The position of the resonances are calculated using 
the most recent experimental data available on 6 Li and are presented in Tab. The 
predicted values of the resonance position are compared with the location we obtained 
experimentally using the following procedure: we prepare 6 Li atoms in the dipole trap 
with the required spin state (mf,m'f) using radiofrequency transfer. We then ramp 
up the magnetic field from to a value B\ slightly higher than the predicted position 
of the Feshbach resonance. The magnetic field is then abruptly decreased to a variable 
value B close to resonance. After a waiting time of 50 ms in the trap, we measure the 
atom number by time of flight imaging. In Fig. 0a, we show the evolution of the atom 
number vs magnetic field in the channel (1/2,-1/2). We observe a sharp decrease of 
the atom number at ~ 186.5 G, close to the predicted value 185 G. The two other 
channels display similar losses in the vicinity of the theoretical position of the Feshbach 
resonances (Tab. 0. Note that in this table the experimental uncertainty is mainly 
due to the magnetic field calibration. 

One of the main issue of the physics of Feshbach resonances is related to the lifetime 
of the molecules, and more generally of the atoms, at resonance. Indeed, one of the 
key element that led to the experiments on the BEC-BCS crossover was the increase 
of the lifetime of molecules composed of fermions close to resonance. To address this 
issue, we have measured the time evolution of the atom number N in the sample at 
the three Feshbach resonances. Since the one-body lifetime is ~ 30 s, much longer 
than the measured decay time (~ 100 ms), we can fit the time evolution using the 
rate equation 
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Table 1: Position and two-body losses parameters of the p-wave Feshbach 
resonances of 6 Li atoms in |/i = 1/2, m^) and I/2 = 1/2, m/ 2 ) . 



^ = -G 2 W-L 3 <n 2 ), (1) 

where n is the atom density and (n a ) — J d?rn a+1 /N (a — 1,2) is calculated from 
the classical Boltzman distribution. 

In contrast to s-wave Feshbach resonances where dipolar losses are forbidden in the 
/ = 1/2 manifold |33| . the losses near a p-wave resonance are found to be dominantly 
2-body in the (1/2,-1/2) and (-1/2,-1/2) channels. The variations of the 2-body loss 
rate with temperature are displayed in Fig. [3]c and show very different behaviors in 
these two channels. This non trivial dependence is actually not the consequence of 
some specific property of the states involved (eg., the quantum statistics) but can be 
recovered using a very simple three-state model. We describe inelastic processes by 
two non interacting open channels, respectively the incoming and decay channels, that 
are coupled to a single p-wave molecular state . This model leads to a very simple 
algebra that can be treated analytically 1341 and yields for the two-body loss-rate at 
a given energy E 

»W = (E-W + W (2) 
Here S = fj,(B — Bp) is the detuning to the Feshbach resonance and K, \i and 7 are 
phenomenological constants depending on the microscopic details of the potential. 
For each channel, these parameters are estimated from our numerical coupled-channel 
calculation (Tab. 0. 

Eqn. 42) shows that, contrarily to s-wave processes that can be described accurately 
by their low energy behavior, p-wave losses are dominated by the resonance peak lo- 
cated at E — S where the losses are maximum. In other word, the so-called "threshold 
laws", that give the low energy scattering behavior, are insufficient to describe a Fes- 
hbach resonance associated with p-wave molecular states (and, more generally, any 
non zero angular momentum molecular state). 

To compare with experimental data, Eq. is averaged over a thermal distribution 
and for 8 > and i5 > 7 we get: 

G 2 ~V^fAy /2 e-* /feBT (3) 



7 V^bT, 

Eqn. |2]is used to fit the data of Fig. |3]b, with B — Bp as the only fitting parameter. 
We get a fairly good agreement if we take B — Bp = 0.04 G (resp. 0.3 G) for the (-1/2,- 
1/2) (resp. (1/2,-1/2)) channel, illustrating the extreme sensitivity of G2 to detuning 
and temperature. This feature was also tested by measuring the variations of G2 with 
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magnetic field at constant temperature. Note that the width of the resonance, as 
given by Eqn. |H] is of the order of k^T / 'fi. At ~ 10 fiK, this yields a width of ~ 0.1 
G, which is comparable with the one observed in Fig. |3]a. 

In the (1/2,1/2) channel, dipolar losses are forbidden and we indeed find that 3-body 
losses are dominant. The dependence of L3 with temperature is very weak (Fig. |3]b) 
and contradicts the low energy theoretical calculation of the temperature dependence 
of three-body loss rate performed in 1351 . Indeed, Wigner threshold law predicts that 
at low energy, L3 should be proportional to T A , with A > 2 for indistinguishable 
fermions. However, as we noticed earlier in the case of two-body losses, the threshold 
law is probably not sufficient to describe losses at resonance due to the existence of 
a resonance peak that might also be present in three-body processes. This suggests 
that 3-body processes must be described by a more refined formalism, analogous to 
the unitary limited treatment of the s-wave elastic collisions 36 . 

Finally, we have checked the production of molecules in (l/2,-l/2)mixture by using 
the scheme presented in 1131 114| . We first generate p-wave molecules by ramping in 
20 ms the magnetic field from 190 G> Bp to B nuc = 185 G < Bp. At this stage, 
we can follow two paths before detection. Path 1 permits to measure the number Ni 
of free atoms: by ramping down in 2 ms the magnetic field from 185 G to 176 G, 
we convert the molecules into deeply bound molecular states that decay rapidly by 
2-body collisions. Path 2 gives access to the total atom number N2 (free atoms + 
atoms bound in p-wave molecules). It consists in ramping up the magnetic field in 

2 ms from B nuc to 202 G> Bf to convert the molecules back into atoms. Since 
the atoms involved in molecular states appear only in pictures taken in path 2, the 
number of molecules in the trap is (N2 — Ni)/2. In practice, both sequences are started 
immediately after reaching B nuc and we average the data of 25 pictures to compensate 
for atom number fluctuations. We then get ATj = 7.1(5) x 10 4 and N 2 = 9.1(7) x 10 4 
which corresponds to a molecule fraction 1 — N1/N2 = 0.2(1). Surprisingly, we failed 
to detect any molecule signal when applying the same method to (1/2,1/2) atoms. 

3 Heteronuclear Feshbach resonances 

So far, the Feshbach resonances used in this paper were involving atoms of the same 
species (namely 6 Li). However, it was recently pointed out that the observation of 
Feshbach resonances between two different atom species could lead to a host of inter- 
esting effects ranging from the observation of supersolid order |37| to the study polar 
molecules |38|. In the case of a mixture of bosons and fermions, these molecules are 
fermions and are expected to be long lived. Indeed, Pauli principle keeps molecules 
far apart and prevents inelastic collisions. Such resonances were observed experimen- 
tally in 6 Li- 23 Na 39 and 40 K- 87 Rb HJ mixtures. In the case of 6 Li- 7 Li in the stable 
Zeeman ground state \ f = 1/2, m/ = 1/2) (g> |/ = l,m/ = 1), the existence and the 
position of heteronuclear Feshbach resonances were predicted in |41| . In that work, the 
6 Li- 7 Li interaction potential was extracted from the data on 6 Li- 6 Li scattering prop- 
erties by mean of a simple mass-scaling. Using coupled channels calculation it was 
found that this system exhibited five Feshbach resonances whose position is displayed 
in Tab. I2l|l2l. 

Experimentally, we probed these Feshbach resonances using a mixture of respec- 
tively Ne, ~ N7 ~ 10 s atoms of 6 Li and 7 Li in the absolute Zeeman ground state 
1 1/2, 1/2} Cg> |1, 1). The gas is cooled at a temperature of 4p,K in the cross dipole 
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B t h (G) B cxp (G) 

218 not seen 

230 226.3(6) G 

251 246.0(8) G 

551 539.9(8) G 

559 548.6(9) G 



Table 2: Predicted and observed positions of the heteronuclear Feshbach resonances 
of a 6 Li- 7 Li gas. First column: resonance position as predicted by coupled channel 
calculation. Secund column: experimental determination of the resonances. The 
experimental error bars are due to the systematic uncertainty on the calibration of 
the magnetic field. The first Feshbach resonance is extremely narrow and could not 
be detected with our experimental resolution. 

trap and we locate the resonance by sweeping down in 1 s the magnetic field from a 
value located about 20 G above the predicted position of the resonance to a variable 
value B. By looking for the value of B at which we start loosing atoms, we were 
able to detect four of the five resonances at a value very close to that predicted by 
theory (Tab. |5J. The discrepancy between the experimental and theoretical values is 
larger in this case than in the case of th p-wave resonance. This is probably due to a 
breakdown of the Born-Oppenheimer approximation that one expects in the case of 
light atoms such as lithium and that forbids the use of the mass scaling 1411 . Note 
also that the missing resonance is predicted to be very narrow (~1 mG). 

As noticed in 1 13| . the width of a Feshbach resonance strongly influences the 
molecule lifetime that can be much higher close to a wide resonance. Since the Fes- 
hbach resonances we found in 1 1/2, 1/2} ® |1, 1) are very narrow (~ 0.1 G) it might 
prove interesting to look for wider Feshbach resonances. In the case of 6 Li- 7 Li, a wide 
resonance is predicted to exist in the stable state 1 3/2, 3/2} ® |1, 1} at ~ 530 G, in 
agreement with our observation of a large atom loss located between ~ 440 G and 
~ 540 G. 
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Figure 1: a,b: Onset of Bose-Einstein condensation in a cloud of 2 x 10 4 6 Li dimers 
at 770 G (a) and of 2 x 10 4 7 Li atoms at 610 G (b) in the same optical trap, (a): 1.2 ms 
expansion profiles along the weak direction x of confinement, (b): 1.4 ms expansion. 
The different sizes of the condensates reflect the large difference in scattering length 
a m = 170 nm for 6 Li dimers and 07 = 0.55 nm for 7 Li. Solid line: Gaussian+Thomas- 
Fermi fit. Dashed line: gaussian component. Condensate fractions are 44% in (a) and 
28% in (b). c,d: BEC-BCS crossover, (c): anisotropy of the cloud, (d): release energy 
across the BEC-BCS crossover region. In (d), the dot-dashed line corresponds to a 
T = ideal Fermi gas. The dashed curve is the release energy from a pure condensate 
in the Thomas- Fermi limit. The solid curve corresponds to a finite temperature mean 
field model with T = 0.5 T<2. Arrow: fc F ja| = 3. 
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Figure 2: Molecular condensate loss rate j3 as a function of the atomic scattering 
length a near the 834 G s-wave Feshbach resonance. The line is a power law fit with 
f3 ~ a" 1 - 9 * - 8 i n agreement with theory, (3 ~ a" 2 55 (T5| . 




Figure 3: (a) Atom number vs. magnetic field £?o,f after a 50 ms wait for atoms 
in the spin mixture (1/2, —1/2) at T ~ 14/^K. The sharp drop close to _Bo ~ 186 G 
over a range ~ 0.5 G is the signature of the p-wave Feshbach resonance predicted by 
theory, (b) (resp. (c)) Variations of 3-body (2-body) loss rates vs temperature at the 
Feshbach resonance, (b): ♦: atoms in the Zeeman ground state |/ = 1/2, m/ = 1/2), 
B ,f ~ 159 G. (c): ■: atoms polarized in \f = 1/2, mj = -1/2), S ,f ~ 215 G. •: 
mixture \f = 1/2, mf = 1/2) + \f = l/2,m f = -1/2), B ,£ ~ 186 G. In both cases, 
the full line is a fit to the data using prediction of Eq. [3] with the magnetic field as 
the only fitting parameter. 
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